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>< Abstract 

We present a new method solving Schwinger-Dyson equations in gauge theories via imaginary 
part analysis. The self-consistency equation for the fermion self-energy can be reduced from 
a numerical unstable eigenvalue problem to a finite dimensional integral which can be solved 
easily. We can therefore derive relations between the dynamically generated fermion mass, 
the coupling, the gauge boson mass and the cutoff. For infinite cutoff the fermion mass 
is independent of the gauge boson mass which should be the relevant scale. Thus for a 
defined fermion mass it is necessary to incorporate a running coupling constant. Within our 
framework this is possible in a well defined way. The results are reasonable and in agreement 
with the renormalization group in the asymptotic regime. 
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1 Introduction 



Almost every theory describing physics beyond the electro-weak scale incorporates some kind 
of dynamical scenario. Although symmetries play an important role in quantum field theory 
and are well understood, most of them are broken and the breaking can only effectively be 
described. Either the Higgs is a composite particle with a rich dynamical background or it is 
replaced by a dynamical model like technicolor. Even if the Higgs is a fundamental particle, 
we need some supersymmetric extension of the standard model to get rid of the quadratic 
divergences. But then the SUSY breaking presumably is caused by some underlying dy- 
namics. Whenever a symmetry is broken a dynamical mechanism is proposed. One reason 
might be the fact that up to now only a running coupling can mediate a symmetry breaking 
scenario with a large hierarchy of scales by blowing up or vanishing at a certain point. 

In contrast our possibilities to deal with dynamical symmetry breaking are restricted. Most of 
the statements stem from simple models (e.g. Nambu-Jona-Lasinio, Gross-Neveu, etc.|l|) or 
from renormalization group arguments 0. However the predictions of dynamical models, e.g. 
gauge models of top condensation [Q], require more effort to find appropriate approximations 
and reliable numerical analyses. 

In this paper we concentrate on scenarios where dynamical symmetry breaking is mediated 
by strongly interacting massive gauge bosons. The corresponding Schwinger-Dyson equa- 
tionsQ are hard to solve within a reasonable approximation. One is confronted with a lot 
of problems. Even in ladder approximation the main parameters like the critical coupling 
and the dynamically generated mass of the fermion can only be calculated with large un- 
certainties. The solutions of the integral equation strongly depend on the behavior at large 
momenta. The eigenvalue equation thus involves numerical instabilities especially for infinite 
cutoff. Furthermore the integral equations are usually solved in the euclidean region. Due 
to the above mentioned problems and the cut in the complex plane it is nearly impossible 
to make statements about the behavior in Minkowski space by analytic continuation. 

One of the key features of dynamical symmetry breaking in asymptotically free field theories 
is that the running coupling defines the scale of chiral symmetry breaking by dimensional 
transmutation. However a full Schwinger-Dyson system with a self-consistent equation for 
the vertex is rather complicated and needs some simplifications. Even if one uses the per- 
turbative renormalization group running for the coupling, the infrared singularity normally 
requires assumptions about the behavior of the coupling in the confinement region. 

^ For an overview see e.g. M. 
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We present here a new method which solves most of the above problems. The integral 
equation is directly treated in Minkowski space, which simplifies the equation, which is then 
solvable by integration. This procedure also enables us to introduce a running coupling in an 
elegant way, which is well defined and does not require any assumptions about the behavior 
in the nonperturbative infrared region. 



2 The Imaginary Part Method 

We consider Schwinger-Dyson equations in the framework of the ladder approximation. The 
self-consistency equations for the vertex function and the vacuum polarization of the gauge 
boson are disregarded. The integral equation for the fermion self-energy is restricted to 
single gauge boson exchange and the above mentioned quantities are substituted by their 
treelevel values. (A running coupling will be incorporated in chapter The resulting 
integral equation is often called the massive gap equation: 




(1) 



-'^ - l/|^(rv).(.)(rv) -'(;^-;3f" . (2) 



Here the includes the gauge part of the coupling. We use the Landau gauge because the 
renormalization constant v4(p^) of the full fermion propagator 

tS-\p) = A{p')^ - Eip') (3) 

is equal to one in that case. With this we get 

y(v'^) = iWT" [ ^^^'^^ ^ m 

' 7 (27r)4p- [S(P)]2(j9-A;)2-M2 ■ ^ > 

The physical solution of this integral equation has to vanish for momenta p'^^oo, otherwise 
the integral is divergent. Therefore the self-energy function cannot be analytical in the 
complete complex plane unless it is trivial. This is also reasonable since for momenta p>M{m 
the self-energy function S(p2) should have an imaginary part because the fermion is then 
able to emit a real gauge boson. We will see this property later in our solutions. For 
simplicity we use the denominator approximation 
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where the S(p^) -function in the denominator is replaced by a fixed mass m. As long as this 
fermion mass is below the boson mass M the S(p^) in the denominator does not induce any 
dynamics and can therefore be replaced[[ 

Up to that point we have used the conventional procedure. This equation is usually in the 
euclidean region where S(p^) is strictly real. In Minkowski space one expects a coupled 
system of real and imaginary part, which seems hard to solve. But we will see that this is 
not the case. We can establish a decoupled integral equation only for the imaginary part, 
first shown in The real part is then calculable by a dispersion relation. 

First we express ^{p^} through the dispersion relation 

_i J 

TX 



+ te 



(6) 



where /(p^) = 3'm[S(p^)] and introduce Feynman parameters. Angular integration results 
then in the one-dimensional integral equation 
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(7) 



with the integral kernel 
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(8) 



where p and I denote the positive square roots of p^ and f^. 

It can now easily be seen that the imaginary part /(p^) starts at M+m: Let I{p^) vanish for 
p^K /i^. The first part of the kernel contributes for p> M+m and the second part for p>M+/i. 
Consider fi<M+m, then only the second part of the kernel contributes for p<M+m. Hence 
/i = M+/i must be true, which is however impossible for M>0. Thus /(p^) vanishes below 
M+m. As mentioned above this is exactly the threshold where a real boson can be emitted. 
One should also introduce a cutoff A to study the dependence on the scale of new physics. 

Using the well known A-function 



X{x,y,z) = \J x"^ + y"^ + z"^ — 2xy — 2xz — 2yz 



(9) 



^For Af = a rich dynamical structure coming from the denominator can appear jsj. 
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and imposing p>M+m we get 

A2 



Hp' 



c 



{p-Mf 



-df 



The first integral is known from tlie dispersion relation: 



(10) 
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The integral equation therefore has the form: 



(11) 
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df 



(12) 



Surprisingly this equation is only formally an integral equation. For M+m<p<2M+m 
the second term in eq.(|l^) vanishes and Hp') is given analytically by the first term. For 
2M+m<p<3M+m Hp') is additionally a simple integral of Hp') o'^^^ already known 
region and so on. Thus Hp') can be expressed by a finite dimensional integral and a finite 
series in C. 

This is an important progress in solving Schwinger-Dyson equations since up to now the 
integral equations were solved as a rather unstable eigenvalue problem. The dynamical mass 
function decreases rather slowly so that an iteration procedure converges badly. On the other 
hand a determination of the eigenvalues of a huge matrix is very extensive and sensitive to 
the stepsize. 

Nevertheless a determination of C for a fixed cutoff A is complicated. One better deter- 
mines A as a function of C, i.e one determines Hp') ^^at scale A where the relation 
(11) is fulfilled. Further integration to the next points where Hp') this normalization 
requirement yields a whole series of possible cutoffs A2, A3, . . .. This is equivalent to different 
eigenvalues C for a fixed cutoff. 



3 Numerical Solutions 

A numerical analysis of eg. ([12|) is shown in fig.|l] and fig.^. We can see that the imaginary 
part of the self-energy function starts at the threshold m + M as we mentioned before. The 
real part has therefore a "resonance like" peak at this point. In fig.0 we multiplied the 
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Figure 1 Real and imaginary part for Figure 2 Asymptotic behavior of the real 
m = 0.5M and C = 0.3. and imaginary part for m = 0.5M, C = 



0.3 and A. = 4 • 10 
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self-energy function with tlie momentum p to compensate tlie 1/p drop of tlie solutions and 
to see the oscillating behavior in the asymptotic region. We took the fourth cutoff, which we 
denoted by A4 to determine the real part of S. The obtained asymptotic behavior is exactly 
what we expect for large momenta as it is shown in appendix ^: 



3?e [S(/) 




S(p') oc-sin A C--ln^ + ^ . (13) 



This demonstrates the numerical stability of the method over a wide range of momenta. 
Additionally we can calculate the cutoff for different couplings C and masses m and M. 
For each value of C one gets a series of cutoffs Aj, for a fixed cutoff one gets a series of 
eigenvalues Ci respectively, which is shown in fig.^. Below the critical coupling Ckrit = 0.25 
dynamical symmetry breaking does not occur, which was analytically shown by Maskawa 
and Nakajima [0. This can now numerically be confirmed by a regression to the point A—*oo. 
This is true for all eigensolutions and astonishingly also for different fermion masses m. But 
this suggests that it is impossible to give a defined fermion mass depending on the coupling 
C for infinite cutoff A, which can be seen in fig.^. The slope of the function m{C) gets bigger 
for rising cutoff. For infinite cutoff it is a vertical line, so the fermion mass is either zero for 
couplings below the critical value or infinite for couplings above. Normally one would expect 
the fermion mass to be dominated by the mass scale of the gauge boson. This is however not 
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Figure 3 C for different cutoffs A with Figure 4 The C- dependence of m for dif- 
m = 0.1 M. ferent cutoffs A. 

the case, the fermion mass is shifted up to the cutoff scale, which corresponds to a hierarchy 
problem.0 The coupling C is large everywhere and a boson mass is not necessary for chiral 
symmetry breaking. So why should the fermion mass be small? The cutoff A is not at all 
removed from the problem although the dynamical fermion mass is finite. We need a running 
coupling in addition to get rid of this problem which is the subject of the next section. 



4 Generalizations — The Running Coupling 

Let us go back one step and have a look at equation (|^) again. This integral equation has a 
much deeper meaning. It can be written in the following form 

= - / (G{p\ M\ m') - G{p\ M\ f )) (14) 



where 

^Thus a serious problem shows up in models where the dynamical symmetry breakdown is caused by a 
broken abelian gauge symmetry. There the coupling is not only constant but rather rises in the UV region. 
So the model has to be embedded in some other theory to remove the Landau pole. The fermion mass is 
then not only pulled up to the U{1) gauge boson mass but to the embedding scale. 
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is the imaginary part of the normal first order perturbation theory. G{p'^,M'^,m'^) can also 
be determined via the optical theorem: 

The inner momenta are therefore strictly on shell. This enables us to introduce a running 
coupling in an elegant way. The full vertex depends on the squares of the momenta, i.e. p^, 
and m^. Since > M^,m^ we can use the mass independent running coupling as an 
ansatz for the vertex. In this case is the only relevant scale which dominates the running 
of the vertex. The same arguments can be applied to G{p'^,M'^,i'^). A running coupling 
can therefore be introduced in a very elegant and well controlled approximation by simply 
replacing C by C(p^).0 This appears to be a very trivial transformation, but the results are 
not as trivial. A very nice feature of our calculation is that we do not have to make any 
assumptions about the coupling below the scale M+m. Although for an asymptotically free 
theory the perturbative coupling develops a pole Air in the infrared, as long as M+m > Air 
we need not care about the pole in the coupling and the change of the running below the 
breaking scale and thus avoid nonperturbative questions in many cases. We have therefore 
a kinematical cut shown as dashed line in fig.| and fig.^ 

To determine the running coupling C(p^) we define the group theory factors Ci, C2 and C3 of 
a vector gauge theory (in brackets the values for the fundamental representation of SU(A^)): 

facdfbcd = Ci5ab (Cl = A^) 

TriT^n) = C26ab (c2 = I) (17) 

ETaTa = Csl (C3 = ^ 

a ^ 

With this definitions we get 

The one-loop running coupling can be written in the form 



Thus we obtain 



9\p') = -^^ with /5o = ^CI-|Ec^ (19) 

'ha 

C(p2) = ^^ with fi: = ^. (20) 
In ^ 



■^In the full system of integral equations only one vertex is substituted by the full vertex to avoid double 
counting of diagrams. However for the considered one loop diagram one has to substitute both vertices since 
the running coupling does not incorporate the dynamical mass function. 
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In fig.|^ the numerical results for different K-values are shown in the limit A— s>oo. The 
dynamically generated fermion mass depends on the ratio M/Ajr. If the boson mass is 
bigger than a certain scale, the fermion remains massless. This corresponds to a critical 
value of the coupling at the breaking scale. If k is smaller than ^ 1 and M < Air the 
imaginary part gets large contributions from the coupling around the pole at Air. Even if 
K goes to zero a dynamical mass m should unphysically be generated with m >M— Air. 
Because of the logarithmic pole structure of the coupling the difference m — Air goes very 
fast to zero so that we cannot calculate m for very small values of k. There one should thus 
replace the pole of the coupling by a more smoothly running coupling which tends to a finite 
value at zero. The m — k dependency is shown in fig.^and ^. For big values of k the fermion 
mass m depends logarithmically on k. 

SU(A^) gauge theories give k > ^ ^ 0.41 for large and are therefore just in the valid 
region. Theories with large fermion representations like QCD with colorsextet quarks or 
technicolor can have even larger fermion masses than Air because of a larger K-value. For 
QCD including only the three lightest quarks u,d,s, which have current masses smaller than 
the infrared pole, we get K = I 0.44. Since M= we obtain a "constituent mass" around 
AqcdS- 

^For the average values of the current quark masses and as ~ 0.118 one gets Aq^^ = 327 MeV. (For the 
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Figure 7 log (x^) for different K-values Figure 8 ^ for small values of k and 
and M = 0. M = 0. For the dashed line see text. 

Note that we do not need any further input parameters to determine this relation between 
Air and the constituent quark mass provided that k is large enough. It is calculated from 
first principles and only depends on the approximation. 

As in the case without running coupling, we can get more than one solution, but now a finite 
number of solutions (See fig.^]). These "excited solutions" are particle excitations related to 
different vacuum states with smaller masses and can therefore only appear in a bubble of 
the corresponding vacuum. It is not clear, whether there is any physical relevance of these 
additional solutions. However if one plays with this consideration, as it is done in P], the 
running coupling provides the possibility to have a finite spectrum, which is controlled by 
M and Air. 

As in the non running case we can compare the asymptotic of our solution, now with the 
asymptotic which is given by the renormalization group. One expects for the irregular and 
regular part of the solution (||T0| and appendix 0) 



^irreg ^ 



and 



S]-gg oc 



(21) 



A fit to the asymptotic part of our solutions reproduces the anomalous dimension within 
a few percent. We only see the regular solution up to scales of the cutoff. This is clear. 



used three-loop beta- function see S) 
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because we did not introduce a current mass up to now. This follows in the next chapter. 



5 Schwinger— Dyson Equation with a Current Mass 



So far we have assumed that the integral 

A2 



/ 



^^^'^ de' (22) 



k^-i^ + IE 

is convergent. In that case we have automatically obtained only regular solutions 

(Inp)** 

^res OC . (,2oj 

This is consistent with the assumption since the integral (|2^) is convergent. If it diverges 
the dispersion relation (^) is not correct. We need then in addition the integral over a half 
circle in the complex plane which compensates the divergence. This infinite constant can be 
replaced by an integral along the real axis as it is done for subtracted dispersion relations 

= " / V. .- 'f' + r / ../'^ , + • (24) 



TT J k"^ — P + ie IT J —u? — P + ie 

— oo — oo 

The last two terms are arranged in such a way that S(— /i^) = m^. Thus we have the 
freedom to fix the mass function S at a certain euclidean]^ scale and therefore to choose a 
relatively large fermion mass at scales far beyond the dynamical symmetry breaking scale. 
This corresponds to having a bare mass, which is the sum of the last two terms in eq.(P^ and 
is, of course, infinite. Putting eq.(0) into eq.(|D leads to the same equation (p!2D. However 
instead of eq. ([Tl|) we have another cutoff-mass relation 



A2 



mvr 



(M+mY ^ 



df + m^TT . (25) 



Using the solution /(p^) for given masses one can calculate the cutoff Aq from eq.(0), which 
is smaller than the cutoff A^^ from eq. (^5|) . Between both cutoffs /(p^) has an irregular 
running which corresponds to the renormalization group running. In the case of a bare mass, 
the cutoff Aq corresponds therefore to the scale where chiral symmetry breaking occurs. 

^ Since the fermion mass is getting complex beyond the threshold it is useful to define it at an euclidean 
scale. 
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6 Conclusion 



In this work we have presented a new method solving Schwinger-Dyson equations in ladder 
approximation. We have established the corresponding integral equation for the imaginary 
part of the dynamical mass function, which is no more an eigenvalue problem but can be 
solved by integration. The procedure is therefore numerically stable. The perfect agreement 
of the shape to the analytic solution in the asymptotic regime is an impressive verification for 
this. Wc have calculated mass - coupling relations for the non-running coupling case, which 
states whether there is a finetuning or not. If the cutoff goes to infinity the fermion mass 
also diverges in the broken phase and thus no defined value is predicted. This is even worse 
than the usual finetuning problem in Nambu-Jona-Lasinio models between the fermion mass 
and the Fermi scale. If the contact interaction represents the exchange of a massive gauge 
boson the Fermi scale is no longer the cutoff of the gauge theory but is connected to the 
gauge boson mass. In the full theory the fermion mass is not shifted to the Fermi scale but 
to the cutoff. The only way out of this dilemma seems to be the introduction of a running 
coupling within an asymptotically free theory. We have shown that in our framework this 
is possible in a well defined way. Now the fermion mass is calculable in dependency of M, 
Air and k and for A— >oo. These are only very fundamental constants of a gauge theory and 
in principle no further assumptions are necessary. Our results are in good agreement with 
the common knowledge about dynamical symmetry breaking, only for very small values of 
K more assumptions seem to be necessary. This should be a point of further investigations. 
The application to many questions in models of electro-weak symmetry breaking like e.g. 
technicolor, topcolor or combined models seems now to be possible with the here developed 
tools. 

Acknowledgments: We would hke to thank Richard Dawid, Dimitris Kominis and Manfred 
Lindner for inspiring and clarifying discussions. 
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Appendix 



A Asymptotic Properties 

The asymptotic properties of the Schwinger-Dyson equation are identical to the massless 
hmit M= since the boson mass M has no influence on the general shape of the solution 
at large momenta. In that limit the integral equation can be transformed to the differential 
equation 




which has the solution [|Tl|] 

S(p^) = m..f,(i + /i^.i-/i^;2;^) . (27) 

Here 2F1 is the hypergeometric function. Its asymptotic can also be derived directly from 
the differential equation (p6|) by neglecting m in the denominator and using the ansatz 
S oc One finds 

52_^ + C' = or ^=^± . (28) 

Only if C is bigger than 1/4 the solution of the differential equation also fulfills the integral 
equation. Otherwise the boundary conditions cannot be fulfilled. Hence 5 is complex. The 
real and imaginary part of the asymptotic has therefore the form 




B Asymptotic Properties for a Running Coupling 

In the case M = our integral equation has the form 



/(/) = C{p^) I mvr^ - / de^^^ I (30) 



which can be transformed to the differential equation: 



^dp^ {^^^'^c\p^)) ~ p^-L • ^^^^ 
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Using the running coupling 



we make the ansatz 



2\-S 



In 



P 



A 



IR, 



(32) 



to describe the asymptotic properties. We consider the leading terms in and In ^ and 
find the relations 6{1 — 6) = and (1 — 25) (1 — 7) = — which give the two solutions 



J(p2) oc 



1 



and /(p^) oc 



M) 

p2 



2 \ K—1 



(33) 



which correspond to the real counterparts (pl|). 
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